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Abstract. 

We present a measurement of the relative branching ratio of the decay — ir^e^iy^ {K^s^) with 
respect to — TT^e^u (7) {Ke3 + KeSj) dccay. The result is based on observation of 19 000 

Kes-y and 5.6 • 10^ K^s decays. The value of the branching ratio is Br{K^r^^, E* > 30 MeV, 6**^ > 
20°)/Br(K°3) = (0.964 ± 0.008l{];[j^^)%. This result agrees with theoretical predictions but is at 
variance with a recently published result. 



1 Introduction 

The study of radiative Ki^ decays can give valuable information on the kaon structure. It allows a 
good test of theories describing hadron interactions and decays, like Chiral Perturbation Theory 
(ChPT). Here we present a study of the radiative Kgs decay. 

There are two distinct photon components in the radiative K^^ decays - inner bremsstrahlung (IB) 
and direct emission. K^^ decays are mainly sensitive to the IB component because of the small 
electron mass rrig. A big contribution to the rate, dominated by the IB amplitude, comes from the 
region of small photon energies E* and angles 6*^ between the charged lepton and the photon, 
with both E* and 0*^ measured in the kaon rest frame. K^^^ amplitude has infrared singularity 
at 0, and a collinear singularity at 6^^ — when rUe = 0. For this measurement and the 

corresponding theoretical evaluation, we exclude these regions by the restriction E* > 30 MeV and 

ei^ > 20°. 

Two different theoretical approaches for evaluation of the branching ratio have been used. Current 
algebra technique together with the Low theorem were applied by Fearing, Fischbach and Smith 
(called FFS hereafter) [1],[2] and by Doncel [3]. ChPT calculations were performed in [4], [5] and 
are being continuously improved [6], [7]. The ratio of the K^^^ to K^^ decay probabilities, applying 
the standard cuts on E* and 9*^, is predicted to be between 0.95 and 0.99%. The amounts of direct 
emission in these various calculations differ, and are roughly estimated to be between 0.1 and 1 % 
of the size of the IB component. 

Two experimental measurements of the K^^^ branching ratio have been published. The NA31 exper- 
iment obtained 5r(i^03^„^; > 30 MeV, 0*^ > 20°)/Br(K03) = (0.934±0.036l[]:[j|)% [8]. The KTeV 
experiment gave a compatible value of the ratio Br{K^^^,E* > 30 MeV, 0*^ > 20°)/Br(Kg3) = 
(0.908 ± 0.008lg:[j}^)% [9]. However, this value does not agree well with theoretical predictions. 



2 Experimental setup 

The NA48 detector was designed for a measurement of direct CP violation in the system. Here 
we use data from a dedicated run in September 1999 where a beam was produced by 450 GeV/c 
protons from the CERN SPS incident on a beryllium target. The decay region is located 120 m 
from the Kl target after three collimators and sweeping magnets. It is contained in an evacuated 
tube, 90 m long, terminated by a thin (3 • 10^^ Xq) kevlar window. 

The detector components relevant for this measurement include the following: 
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The magnetic spectrometer is designed to measure the momentum of charged particles with 
high precision. The momentum resolution is given by 



= (0.48 e 0.009 • p) % 



(1) 



P 



where p is in GeV/c. The spectrometer consists of four drift chambers (DCH), each with 8 planes 
of sense wires oriented along the projections x,u,y,v, each one rotated by 45 degrees with respect to 
the previous one. The spatial resolution achieved per projection is 100 /im and the time resolution 
is 0.7 ns. The volume between the chambers is filled with helium, near atmospheric pressure. The 
spectrometer magnet is a dipole with a field integral of 0.85 Tm and is placed after the first two 
chambers. The distance between the first and last chamber is 21.8 m. 

The hodoscope is placed downstream of the last drift chamber. It consists of two planes of scin- 
tillators segmented in horizontal and vertical strips and arranged in four quadrants. The signals 
are used for a fast coincidence of two charged particles in the trigger. The time resolution from the 
hodoscope is 200 ps per track. 

The electromagnetic calorimeter (LKr) is a quasi-homogeneous calorimeter based on liquid 
krypton, with tower read out. The 13212 read-out cells have cross sections of 2 x 2 cm^. The 
electrodes extend from the front to the back of the detector in a small angle accordion geometry. 
The LKr calorimeter measures the energies of the and 7 quanta by gathering the ionization 
from their electromagnetic showers. The energy resolution is : 



where E is in GeV, and the time resolution for showers with energy between 3 GeV and 100 GeV 
is 500 ps. 

The muon veto system (MUV) consists of three planes of scintillator counters, shielded by iron 
walls of 80 cm thickness. It is used to reduce the Kl — Tr^fi^u background. 

Charged decays were triggered with a two-level trigger system. The trigger requirements were two 
charged particles in the scintillator hodoscope or in the drift chambers coming from the vertex in 
the decay region. 

A more detailed description of the NA48 setup can be found elsewhere [10] . 



3 Analysis 

3.1 Event selection 

The data sample consisted of about 2 TB of data from 100 million triggers, with approximately 
equal amounts recorded with alternating spectrometer magnet polarities. These data are the same 
which were used for the measurement of the branching ratio [11]. The following selection criteria 
were applied to the reconstructed data to identify Ke3 decays and to reject background, keeping in 
mind the main backgrounds to Kgs, which are Kl — tt^ij^v {Kf^^) and Kl — 7r"'"7r~7r'^ {K^t^): 




(2) 
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- Each event was required to contain exactly two tracks, of opposite charge, and a reconstructed 
vertex in the decay region. To form a vertex, the closest distance of approach between these tracks 
had to be less than 3 cm. The decay region was defined by requirements that the vertex had to be 
between 6 and 34 m from the end of the last collimator and that the transverse distance between 
the vertex and the beam axis had to be less than 2 cm. These cuts were passed by 35 million events. 

- The time difference between the tracks was required to be less than 6 ns. To reject muons, only 
events with both tracks inside the detector acceptance and without in-time hits in the MUV system 
were used. For the same reason only particles with a momentum larger than 10 GcV were accepted. 
In order to allow a clear separation of pion and electron showers, we required the distance between 
the entry points of the two tracks at the front face of the LKr Calorimeter to be larger than 25 cm. 
As a result 14 million events remained. 

- For the identification of electrons and pions, we used the ratio of the measured cluster energy, 
E, in the LKr calorimeter associated to a track to the momentum, p, of this track as measured in 
the magnetic spectrometer. The ratio E/p for a sample of 75 000 pion tracks, selected by requiring 
the other track of a 2-track event to be an electron with E/p > 1.02, is shown in fig. 1. As a 
cross-check pion samples from K2t^ and i^3,r decays were selected giving similar results. Also shown 
in the figure is the distribution for 450 000 electron tracks which are selected from 2-track events 
where the other track is a pion, with 0.4 < E/p < 0.6. For the selection of i^es events, we require 
one track to have 0.93 < E/p < 1.10 (electron) and the other track to have E/p < 0.90 (pion). 11.7 
million events were accepted. 
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Fig. 1. Distribution of the ratio of the shower energy E reconstructed by the LKr and the momentum 
p reconstructed by the spectrometer, for pions (dotted) and electrons (line) from K^s events (see 
text). 

- In order to reduce background from decays, we required the quantity 

{m%; - m\_ - mlof - 4(m^_m^o -|- m|^p^) 



4(Pj_ + m"^ 



(3) 



+- 



to be less than — 0.004(GeV/c)^. In the equation above, p± is the transverse momentum of the two 
track system (assumed to consist of two charged pions) relative to the flight direction and m+_ 
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is the invariant mass of the charged system. The variable Pq^ is positively defined if the charged 
particles are pions from the decay K^.,^ and its distribution has maximum at zero. The cut removes 
(98.94±0.03)% of K^^ decays and (1.03±0.02)% of K^z decays as estimated with the Monte Carlo 
simulation (sect. 3.3). After this cut, we were left with 11.4 million i^eS candidate events. 

The neutrino momentum in K^,^ decays is not known and the kinematic reconstruction of the kaon 
momentum from the measured track momenta leads to a two-fold ambiguity in the reconstructed 
kaon momentum. The solution with larger energy we call "first solution" . In order to measure the 
kaon momentum spectrum, we selected events in which both solutions for the kaon momentum lie 
in the same bin of width 8 GeV. These 4 • 10^ events we call "diagonal events". 

The last selection criterion was the requirement that each of the two solutions for the kaon energy 
had to be in the energy range (60,180) GeV. As a result of this selection, 5.6-10^ fully reconstructed 
events were selected from the total sample. These selected events include radiative K^^ events. 

For the selection of i^e37 events, the following additional requirements were made: 

The distance between the 7 cluster and the pion track in LKr had to be larger than 55 cm in 
order to allow a clear separation of the 7 cluster from pion clusters. As is shown in fig. 2 the 
hadron showers can extend over lateral distances of up to 60 cm from the track entry point in LKr. 
After the requirements for E* > 30 MeV and 9* > 20° (for both solutions of the kaon energy), 



events/2cm 




1 r 



10 20 30 40 50 60 
distance, cm 

Fig. 2. Transverse distance between the pion entry point in LKr and the position of a cluster induced 
by pion interactions with matter; the pions here are selected from — tt+tt" decays where the 
entry points of the two tracks in the LKr calorimeter are at least 80 cm from each other; clusters 
have a minimum energy of 4 GeV 

22 100 events survived. To distinguish the 7 from the electron cluster we required the transverse 
distance between the 7 cluster candidate and the electron track in LKr to be greater than 6 cm. 
The electromagnetic transverse rms shower width in LKr is 2.2 cm. An event was rejected if the 
7 cluster candidate was less than 16 cm away from the beam axis, because of the beam hole in 
the LKr calorimeter. We also rejected events with a 7 cluster candidate with energy below 4 GeV 
because the energy resolution deteriorates below this threshold. Finally an event was rejected if 



7 



the 7 was not in-time (more than 6 ns time difference) with the associated cluster (s). These cuts 
provided a sample of 19 117 Ke^^ candidates. 

3.2 Backgrounds 

The amount of background was evaluated using a Monte Carlo simulation for other kaon decays. 

The background to ii'e37 events is small and comes from three sources - K^t^ and Kl ir^TT^e^u 
{K^a) decays as well as K^i decays with an accidental photon. The K^j^ background was reduced 
by the cut on the variable and the electron identification through the E/p > 0.93 condition. 
Variations of these cuts have a negligible effect, since the probability to misidentify a pion for an 
electron is only 0.57% from fig. 1, and the Pq^ distribution is well reproduced by the MC simulation. 
The estimated number of background events was 40112 events. 

The Ke4 background was evaluated to be 80 it 40 events from the measured branching ratio and 
the calculated acceptance for these decays. 

The contamination from K^^ decays with an accidental photon was estimated using the distribution 

of the time difference between the 7 cluster candidate and the (average) time of the other cluster (s). 
The number of events in the two control regions (—25,-10) ns and (10,25) ns were extrapolated 
to the signal region (—6, 6) ns. The final number for this source of background was estimated to be 
20I20 events, assuming a flat distribution. 

All backgrounds to Kg^^ add up to 140 ± 82 events or 0.7% of the total KgSj sample of 19117 
events. 

The main background to the normalization channel K^.^ arises from K^j^ and K^^-^ decays. The 
estimations were made as in the case of i^e37- background decays together gave a i^es 

signature in less than 9 • 10~^ of the cases (< 500 events). This percentage is negligible compared 
to background sources in Kg^j decay. 

3.3 Monte Carlo Simulation 

In order to calculate the geometrical and kinematical acceptance of the NA48 detector, a GEANT- 
based simulation was employed [10]. The kaon momentum spectrum from sect. 3.1 was implemented 
into the MC code. The radiative corrections (virtual and real) were taken into account by modifying 
the PHOTOS [12] program package in such a way as to reproduce the experimental data. This was 
achieved by weighting the angular distribution in the centre-of-mass frame such as to fit the 
experimental data (model independent analysis). With this procedure the MC and experimental 
data showed good agreement. As an example, the distributions of the neutrino energy, 7 energy and 
61^ (first solutions) in the centre-of-mass frame are presented in figures 3, 4, and 5 respectively. The 
upper plots of the figures show the experimental data distributions and the lower show the ratio 
of the data and the MC spectra, normalized to unity. The plots represent data with the negative 
magnet polarity and after the K^^^ selection. 

The MC data were treated exactly in the same way as the experimental data and were used for 
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Fig. 3. Reconstructed neutrino energy in the Center of Mass System; upper part - experimental 
data distribution, lower part - normalized to unity ratio of DATA/MC linearly fitted 
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Fig. 4. First solution for E*; upper part - experimental data distribution, lower part - normalized 
to unity ratio of DATA/MC linearly fitted 

acceptance calculations. The acceptance for Kg^^ is e(i^e3-y) = (6.08 ± 0.03)% as compared to the 
Ke3 acceptance eiK^s) = (17.28 ± 0.01)%. 



3.4 Reconstruction and analysis technique 



We used the "diagonal events" to measure the kaon momentum spectrum from K^s decays. However, 
as this reduces the data sample significantly, for the analysis of the branching ratio the problem 
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Fig. 5. First solution for 9*^\ upper part - experimental data distribution, lower part - normalized 
to unity ratio of DATA/MC linearly fitted 

was dealt with in another way. In the K^^ selection it was required that both solutions were in the 
range (60,180) GeV. Further in the K^.^^ selection events were rejected if (at least) one of the two 
solutions for E* was less than 30 MeV or (at least) one of the two solutions for ^*,, was less than 
20". The same procedure was used for selecting MC events when calculating the acceptance. 



An important issue are radiative corrections. Only the inclusive rate (i^e37 P^^^ ^'^y number of 
radiative photons) is finite and calculable. In our selection we have required only one hard 7 
satisfying E* > 30 MeV and 91^ > 20°. In this way in the final selection events with one "hard" 
7 and any number soft photons are included. Events with two or more hard photons are rejected. 
This loss has to be taken into account by MC in the calculation of the corresponding acceptance. In 
order to check the MC we have compared the number of 7 clusters in the LKr calorimeter predicted 
by the MC with the one in the experimental data. A slight difference has been observed leading to 
a small correction of 0.05% to the branching ratio. We take this into account by a correction factor 
Cm = 0.9995 to the branching ratio. Additionally we have reanalysed our data, requiring at least 
one hard photon i.e. accepting any number of photons. This is the inclusive rate which is finite and 
can be calculated. The result for R agreed within 0.2% with the analysis requiring exactly one hard 
photon. 

The trigger efficiency was measured to be (98.1 ± 0.1)% for K^s decays and (98.1 ± 0.6)% for K^s^ 
decays. 

On the basis of 19117 Kcs-^, candidates with an estimated background of 140 it 82 events and 
5.594 million Kg^ events (including additional photons) after background subtraction, and using 
the calculated acceptances, the branching ratio was computed from the relation: 



The result from 9361 KeSj events and 2.728 million Kg^ events for positive magnet polarity was 
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Fig. 6. Kaon momentum distribution obtained from K^z (line), K27r (open squares) and K-^t^ (circles) 
decays. Arbitrary units on Y-axis 

R = (0.953 ±0.010)% and from 9616 Kgs^ events and 2.866 million K^s events for negative polarity, 
R = (0.975 ±0.010)%, where the errors are statistical. We now turn to the systematic uncertainties. 



3.5 Systematic uncertainties 

Our investigation of possible systematic errors showed that the biggest uncertainty comes from 
the kaon momentum spectrum. In order to determine the influence of this factor we reconstructed 
the experimental kaon momentum distribution from K —>■ tt+tt" and K — > 7r"'"7r~7r*^ decays and 
implemented them in the MC simulation. The shape of the spectrum for the three decays is shown 
in fig. 6. The systematic error from the momentum spectrum was estimated by taking the 3 different 
momentum spectra and calculating the effect of this variation on the acceptance ratio of K^z, and 
Kf.^^. It resulted in an relative uncertainty of ( I3) • 10~^. 

The stability of the result upon the various cuts used in the i^e37 selection was also investigated. 
The cuts were varied in between values which rejected no more than 10% of the events. The biggest 
fluctuations in the branching ratio estimation were taken as systematic errors, and all the errors 
were added in quadrature with a relative result of ±5 ■ 10~^. 

Uncertainties in accidental photon events and in other background contributions are dominated by 
statistics and are not amongst the largest of the systematic errors (( 1^) • 10~^ and ( I3) • 10~^ 
correspondingly) . The influence of the K^z selection cuts to the final result was estimated as in the 
case of ^TeST selection cuts. The quadratic addition of all these relative errors from variations of 
individual selection cuts yielded an inclusive relative error of ±5- 10~^. The value of the form- factor 
A+ in the K^.^ decay was varied between 0.019 and 0.029. The largest fluctuation was taken as a 
relative systematic error - ±1 • 10~^. 

Our estimate of the systematic errors is summarized in Table 1. 
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Source 


lAK/ K 


Kl spectrum 


+6 ir,-3 


i^e37 selection 


1 FT 1 n — 3 


7 accidentals 


+2 1 n-3 


Background uncertainties 


-3 ■ J^O 


ife3 selection 


±5 • 10-3 


Form-factor uncertainties 


±1 • 10-3 


TOTAL 


t^l ■ 10-3 



Table 1 

Relative systematic uncertainties to the Branching ratio 
4 Results and conclusion 



The results are based on 18977 K^^-y and 5.594 • 10^ K^z events. We obtain the following value for 
the branching ratio including the systematic error: 

R = (0.964 ± 0.0081°:°^^)% = {OMAt^omiW^ (5) 



Figure 7 shows this branching ratio compared to theoretical and experimental results. The authors 
of ref. [7] have undertaken a serious effort to estimate the theoretical uncertainties in R, while for 
the earlier theoretical values, this error is not known. These authors obtain R = (0.96 it 0.01)%. It 
appears that our experimental result agrees well with the theoretical calculations [2], [3], including 
the most recent one [7]. However our result is at variance with a recent experiment with similar 
statistical sensitivity [9]. Our measurement, with a 1.5% precision, therefore confirms the validity 
of calculations based on chiral perturbation theory. 



NA48 

KTEV 
NA31 



m Doncel, FFS 

ChPT 0(p6) 

0.875 0.9 0.925 0.95 0.975 1 1.025 1.05 1.075 

BR(K^^,20"',30MeV)/BR(K^ "1""^ 

Fig. 7. Theoretical and experimental results for the radiative K^s branching ratio. The two lower 
entries in the plot are theoretical results. 
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